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Abstract 

We present the 5-dimensional cosmological solutions in the Randall-Sundrum warped 
compactification scenario, using the Goldberger-Wise mechanism to stabilize the size 
of the extra dimension. Matter on the Planck and TeV branes is treated perturbatively, 
to first order. The back-reaction of the scalar field on the metric is taken into account. 
We identify the appropriate gauge-invariant degrees of freedom, and show that the 
perturbations in the bulk scalar can be gauged away. We confirm previous, less exact 
computations of the shift in the radius of the extra dimension induced by matter. We 
point out that the physical mass scales on the TeV brane may have changed signif- 
icantly since the electroweak epoch due to cosmological expansion, independently of 
the details of radius stabilization. 



1. Introduction. The Randall- Sundrum (RS) idea []1| for explaining the weak-scale 
hierarchy problem has garnered much attention from both the phenomenology and string- 
theory communities, providing a link between the two which is often absent. RS is a simple 
and elegant way of generating the TeV scale which characterizes the standard model from 
a set of fundamental scales which are of order the Planck mass (Mp). All that is needed is 
that the distance between a hidden and a visible sector brane be approximately b = 35/Mp 
in a compact extra dimension, y G [0, 1]. The warping of space in this extra dimension, by 
a factor e'^^^, translates the moderately large interbrane separation into the large hierarchy 
needed to explain the ratio TeV/Mp. 

However the RS idea as originally proposed was incomplete due to the lack of any mecha- 
nism for stabilizing the brane separation, b. This was a modulus, corresponding to a massless 
particle, the radion, which would be ruled out because of its modification of gravity: the 
attractive force mediated by the radion would effectively increase Newton's constant at large 
distance scales. An attractive model for giving the radion a potential energy was proposed 
by Goldberger and Wise (GW) P|; they introduced a bulk scalar field with different VEV's, 
Vq and Vi, on the two branes. If the mass m of the scalar is small compared to the scale k 
which appears in the warp factor e~^^y , then it is possible to obtain the desired interbrane 
separation. One finds the relation e~^^ = (vi/vq)^^ . 

An important benefit of stabilizing the radion is that cosmology is governed by the usual 
Friedmann equations, up to small corrections of order p/(TeV)'^ [H. Even with stabilization, 
there may be a problem with reaching a false minimum of the GW radion potential [Q, 
but without stabilization, there is a worse problem: an unnatural tuning of the energy 
densities on the two branes is required for getting solutions where the extra dimension is 
static p, |], a result which can be derived using the (5,5) component of the Einstein equation 
Gmn = i^'^Tmn- Howevcr when there is a nontrivial potential for the radius, V{b), the (5,5) 
equation serves only to determine the shift 5b in the radius due to the expansion, and there 
is no longer any constraint on the matter on the branes. Although this point is now well 



appreciated 131-1131' previously been explicitly demonstrated by solving the full 

5-dimensional field equations using a concrete stabilization mechanism. Indeed, it has been 
claimed recently that such solutions are not possible with an arbitrary equation of state for 



the matter on the branes [12|-[13|, and also that the rate of expansion does not reproduce 



normal cosmology on the negative tension brane despite stabilization |]T4|. Our purpose is 
to present the complete solutions, to leading order in an expansion in the energy densities 
on the branes, thus refuting these claims. 

2. Preliminaries. The action for 5-D gravity coupled to the stabilizing scalar field $ 
and matter on the branes (located at y = and |/ = 1, respectively) is 

S = j d''x^[-^R-K + \d,^d^'^-V[^)) 

+ j ci^a;v^(/:^,o-Fo(<f))|j,=o + / d^x^{C^^^-V^{^))\y=^, (1) 

where is related to the 5-D Planck scale M by = 1/(M^). The negative bulk cosmolog- 
ical constant needed for the RS solution is parametrized as A = — 6A;^/k^ and the scalar field 
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potential is that of a free field, V{^) = im^$^. The brane potentials Vq and Vi can have 
any form that will insure nontrivial VEV's for the scalar field at the branes, for example 
Vi{^) = Xi{^^ - vfY 0. In ref. IQ we pointed out that the choice V^,($) = mi(<l> - ViY is 
advantageous from the point of view of analytic calculability (see also [p!5|] ). 
We will take the metric to have the form 



ds"^ = n'^jt, y)dt^ — a'^{t, y) dxf — y)dy'^ 

i 

= e-2^(*'^)rft2 _ «^(i)2e-2A{t,,) j2 dxl - bit, y fdy\ 



(2) 



where a perturbative expansion in the energy densities of the branes will be made around 
the static solution: 



N{t,y) = A^{y)+5N{t,y)- 
h{t,y) = bo + 6b{t,y); 



A{t,y) = Ao{y)+6A{t,y) 
<^{t,y) = <i>o{y)+S<l>{t,y). 



(3) 



The perturbations are taken to be linear in the energy densities and p of matter on the 
Planck and TeV branes, located at y = and y = I, respectively. 

This ansatz is to be substituted into the Einstein equations, Gmn = '^'^Tmn, and the scalar 
field equation 
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n 
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dt -6a^<l> ) -dyi -a^n$' + ba-^n [V + V^Siby) + Vl6{b{y - 1))] = 0. 



(4) 



Here and in the following, primes on functions of y denote while primes on potentials of 
$ will mean The nonvanishing components of the Einstein tensor are 
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and the stress energy tensor is T„ 
On the branes, is given by 



W(V^($) + A) + 5™$5n$ - \d^^di^g^^n in the bulk. 



Tl = S{by)dmg{Vo + p.,Vo-p.,Vo-p.,Vo-p.,0) 
+ S{b{y - 1)) diag(ri + p,Vi- p, Vi -p,Vi- p, 0) 



(6) 
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At zeroth order in the perturbations, the equations of motion can be written as 



$[; = AA',^', + mX%, (7) 



and the solutions are approximately 

%{y) = Voe-^"">y; Aoiy) = kh,y + ^vl{e-^^''^y - 1) 
where we have normalized Aq{Q) = 0, and introduced 

2 

2 ^ — 




(9) 



The above approximation is good in the limit e -C 1, which is the same regime in which the 
Goldberger-Wise mechanism naturally gives a large hierarchy without fine-tuning the scalar 
field VEV's on the branes: e~*^*° = {vi/voY^'^. For small e, the GW solution coincides with 



an exact solution of the coupled equations that was presented in ref. [|15 . 

3. Perturbation Equations. We can now write the equations for the perturbations of 
the metric, SA, 6N, 6b, and the scalar field, 5$. The equations take a simpler form when 
expressed in terms of the following combinations: 

^ = 6A' - A'^ j- - — T = 6N' - 6A' (10) 
Oq 3 

Further simplification comes from realizing that the perturbations will have the form, for 
example, \E' = p*(t)fi'o(z/) + pif)9i{.y)i so that their time derivatives are proportional to p 
and p*. Below we will confirm that p = —3H{p + p), where H ~ ^/p, ^/p^ is the Hubble 
parameter. Therefore time derivatives of the perturbations are higher order in p and p* than 
are y derivatives, and can be neglected at leading order (except in the (05) Einstein equation, 
where p^^^ is the leading order). Using this approximation, we can write the combinations 
(00), (00) — (ii), (05) and (55) of the Einstein equations as 

4A[,vi/-vi>' = (11) 

-AA',T + r = 2((^)' -^]by^" (12) 

-^T + ^ = (13) 
ao 

A',{4^ + T) + 'jU6^-%6^' + ^',^^^] = [{^)' + ^)by^° (14) 
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In addition, there is the scalar field equation, 

= (4V1/ + T)$;, + + blV"{%)^ 5$ + 4A^5$' + {2blV'{<^o) + 4A'o<|.;) ^ + <|.;^ 

(15) 

Assuming Z2 symmetry (all functions symmetric under y —y), the boundary conditions 
implied by the delta function sources at the branes are 

2 2 

*(t,0) = +^&oP*W; ^(t,l) = -^6op(t) (16) 

o o 
2 2 

T(t,0) = -y&o(p*+P*)(t); T(t,l) = +y6o(p+p)(t) (17) 

5$'(t,2/„) = ^MiL^$[,(t,y„) + (-l)"('|W:(<|.o(t,y^^ (18) 



where in ([TsD n = 0, 1, ?/o = and yi = 1. 

4. Solutions. Naively, it would appear that we have five equations for four unknown 
perturbations, but of course since gravity is a gauge theory, this is not the case. First, we 
have the relation ^[Eq. |11|] +^[Eq. |12| = [Eq. |13|. Furthermore, the (55) Einstein equation 
and the scalar equation can be shown to be equivalent, using (00), (ii) and the zeroth order 
relations (0): [Eq. {^]' - 4:A'oX[Eq. (0)] = $ox[Eq. (|15|)]. So our system is actually 
underdetermined because of unfixed gauge degrees of freedom. To see this more directly, 
consider an infinitesimal diffeomorphism which leaves the coordinate positions of the branes 
unchanged: y = y + f{y), where /(O) = /(I) = 0. The metric and scalar perturbations 
transform as 

6 A 6A + A'J; 6N^6N + A'J 

6b 6b + bof; 5$ ^ 5$ + (19) 

If desired, one can form the gauge invariant combinations 

6h k"^ 6h 
6A' -A^-^-^^'^6^- 6N'-6A- ^16^ - %6<^' + <^f-^ (20) 
Oo 3 60 

the first two are precisely our variables \E' and T and the last one appears in (55) equation. 
In terms of these gauge invariant variables, the system of equations closes. 

It is now easy to verify the following solution from the (00) and (00)- (ii) equations, i.e., 
eqs. ( p!l| - |T^ ). Denoting the warp factor Vt = e"^"'-^-', we find 

^ = Q^fX^^ e'^^^y^ {F{y){n'p + p.) - {Q'p + Q'p^)) (21) 
T = ^^fX^) "'^"^'^ {-Fiym\p + P)+P*+P*) + i^\p + p) + ^'ip*+P*))) (22) 
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where 



F(y) = 1 - (1 - JO^ ^ ^ ^-2kboy 



and the Friedmann equations are 

'ao^ _ SttG 
.aoJ 3 
■ ao \ " do 

^CLqJ CLq 



(23) 

(24) 
(25) 
(26) 



The approximations in eqs. ( ^3]) and (|26D hold when the back reaction of the scalar field on 
the metric can be neglected. 

In the Friedmann equations (|2^p5D , we note that p is the bare value of the energy den- 
sity on the TeV brane, naturally of order M^, while fl^p is the physically observable value, 
of order (TeV)"^. Since p* has no such suppression, it seems highly unlikely that p* should 
be nonzero today; otherwise it would tend to vastly dominate the present expansion of the 
universe. We also point out that these equations are consistent only if energy is separately 
conserved on each brane: p + 3H{p + p) = and p* + 3if(p* + p*) = 0. This can be derived 
directly by considering the (05) Einstein equation, evaluated at either of the branes. The 
equations of state on the two branes are completely independent; there is no relation between 
p/p and p*/p*. 

5. Stiff potential limit. The above solutions are quite general, but they are not 
complete because we have not yet solved for the scalar field perturbation, 5$. This would 
generically be intractable, but there is a special case in which things simplify, namely, when 
the brane potentials Vi{^) become stiff. In this case, the boundary condition for the scalar 
fiuctuation becomes 5$ = at either brane. There is no information about the derivative 
5$' in this case; although 5$ — 0, at the same V"{^) — oo in such a way that the product 
(5$^"($) remains finite, and eq. (|1^) is automatically satisfied. 

Notice that the shift in 5$, eq. (|19|) , respects the boundary conditions on 5$. Moreover, 
$0 is always nonzero for our solution. It is therefore always possible, given some solution 
(5$ which vanishes at the branes, to choose an / such that 5$ becomes zero. This is a 
convenient choice of gauge because it simplifies the equations of motion, and we will make it 
for the remainder of this letter .Q Thus far we have satisfied the (00), (ii) and (05) Einstein 
equations. As noted above, eqs. (|T^) and (|T^) are equivalent, so either one just determines 
the shift in the radius. Using the former, and defining 



Giy) 



lp2Ao(|/) 



Jo J Jo 



-2A, 



°dy 



1-^2 ' 



(27) 



^The above argument is strictly true only for difFeomorphisms which are constant in time, while for our 
problem we need /(t, y) ~ p{t), p*{t). However, the time variation of such an / is of higher order in p and 
so we can neglect it to leading order in the perturbations. 
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we find that 

6b 



12 



bo 2%^ 



n\p - 3p)G + (p, - 3p,){G - A'qC 



n\p - 3p) + n\p, - 3p,) ; (2 



the last expression is found by approximating Aq = kboy everywhere, which means neglecting 
the back reaction. Using the zeroth order solution (|]) for $0, and integrating over y, we can 
obtain the shift in the size of the extra dimension, 

'^j^.^mp^MiiliPi^^ (29) 



We can compare this to the result of ref. |]l6l by using their result for the radion mass, 
mf. ^ {A/3)K^{evokfn^+'^% and the relation kn^ ^ l/M^. Then 

IpSbdy ^ [n\p-3p) + n^ip,-3p.)] 

bo Gkbom^^Mp^ ^ ' 



which agrees with ref. [Q, except for small corrections of order (1 + f2^).0 As is well known, 
the shift in the radion vanishes when the universe is radiation dominated, because the radion 
couples to the trace of the stress energy tensor, which vanishes if the matter is conformally 
invariant. 

6. Implications. Above we focused on the shift in the size of the extra dimension due 
to cosmological expansion, but the more experimentally relevant quantity is the shift in the 
lapse function, n(t, 1), evaluated on the TeV brane. As emphasized in ref. the change 



in n{t^ 1) between the present and the past determines how much physical energy scales on 
our brane, like the weak scale, Mw-, have evolved. The time dependence of Mw is given by 
Mw{t)/Mw{to) = e-'5^(*'i)+'5^(*0'i)+'5^(*'0)-'5^(*0'0). In terms of the variables of the previous 
section, 6N' = ^ + T + A'o6b/bo. We find that 

£ SN'it, y)dy - ^ ((2p + 3p) - e'"'i2p^ + 3p.)) + £ A'^^^^dy (31) 

Interestingly, the new non-56 contribution is present even during radiation domination, is 
parametrically smaller than the radion part only in the matter dominated era, and then only 
if the back reaction is small (e ^ 1). If MpQ ~ 1 TeV, the shift in the energy scale since 
nucleosynthesis is negligible, and this is the only cosmological constraint on 6N. However, 
near the weak scale, the effect could be more interesting. To first order in the physical energy 
density pp = fl'^p, at early times, 

M„.W.M.w(l-^^), (32) 



^Our correction factor is (1 - -fl^), while that of ref. hj is (1 - f^^). 
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assuming that p* = = 0. It is conceivable that k ~ Mp/30 — in fact, the RS model requires 
k < Mp for consistency, so that higher dimension operators in the gravitational part of the 
action do not become important. With these parameters and assuming (7* ~ 100 relativistic 
degrees of freedom and VLMp = 1 TeV, the correction to Mw becomes of order unity at 
a temperature of 130 GeV. Thus the temporal variation in fundamental mass scales might 
have some relevance for the electroweak phase transition and baryogenesis. 

We thank C. Csaki, M. Graesser and G. Kribs for helpful discussions. JC thanks Nordita 
for its hospitality while this work was being finished. 
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